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Ground state and finite temperature properties of a system of coupled frustrated and/or 
dimerized spin-1/2 chains modeling e.g. the CuGeOa compound are reviewed. Special empha- 
sis is put on the investigation of the role of impurity doping. A chain-mean field computation 
combining exact diagonalisations of the chain hamiltonians together with a mean field treat- 
ment of the weak interchain couplings is performed in order to map the microscopic model onto 
a low-energy effective model. The latter describes a 2-dimensional system of effective spin-1/2 
local moments interacting by spacially anisotropic long range spin exchange interactions. An 
extensive study of this effective model is performed by Stocastic Series Expansion Quantum 
Monte Carlo for a wide range of temperatures and impurity concentrations. Interesting scal- 
ing behaviors of the uniform and staggered spin susceptibilities (above a small ordering Neel 
temperature due to a residual 3D coupling) can be interpreted in terms of the formation of 
large clusters of correlated spins carrying a finite magnetization. Such results are reproduced 
satisfactorily by a new Real Space RG enabling to deal with long range interactions in two- 
dimensions. 
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1. Introduction 

Low dimensional gapped quantum magnets have at- 
tracted a lot of interest in condensed matter physics for 
many years. The possibility of doping such systems has 
lead to an extremely rich emerging field. For instance, 
the discovery of the first non-organic spin-Peierls (SP) 
compound CuGcOs 1 and its doping with static non- 
magnetic impurities realized by direct substitution of a 
small fraction of copper atoms by zinc 2 or magnesium 3 
atoms offered a new challenge for the theorist. Replac- 
ing a spin-i in a spontaneously dimerized spin chain by 
a non-magnetic impurity, described as an inert site, re- 
leases a free spin-i, named a soliton, which docs not 
bind to the dopant. 4 The physical picture is completely 
different when a static bond dimcrisation exists and pro- 
duces an attractive potential between the soliton and the 
dopant 4,5 and consequently leads, under doping, to the 
formation of local magnetic moments 4, 6 as well as a rapid 
suppression of the spin gap. 7 However, a coupling to a 
purely one-dimensional (ID) adiabatic lattice 8 does not 
produce confinement in contrast to more realistic models 
including an elastic inter-chain coupling (to mimic 2D or 
3D lattices). 8,9 

In the following, we are going to present through a 
generic model of quasi one-dimensional SP system the 
impurity induced local moment formation as well as the 
tendency towards antiferromagnetic (AF) ordering in the 
thermodynamic limit observed for any non-zero dopant 
concentration. An interesting low temperature scaling 
behavior will also be presented and interpreted in term of 
the formation of large clusters of correlated spins carry- 
ing a finite magnetization. This short review is the con- 
sequence of several previous works 10-12 involving many 
numerical tools: Lanczos Exact Diagonalization (ED) as- 
sociated to chain-mean field theory, Stocastic Series Ex- 



pansion (SSE) Quantum Monte Carlo (QMC) methods 
including long-range interactions, numerical Real Space 
Renormalization Group (RSRG) technique. 

The rest of the paper is organized as follows. In the fol- 
lowing section, we introduce a generic microscopic model 
describing weakly coupled frustrated spin-^ chains, and 
compute its phase diagram without impurity. Then in 
section 3 we concentrate on impurity efects, through the 
confinement mechanism, responsible of the impurity in- 
duced local moment formation. We also study the ef- 
fective interaction between impurities which enable us 
to construct an effective two-dimensional (2D) diluted 
model. The section 4 is thus devoted to the impurity 
induced AF ordering which appear in such an effective 
model. Finally, in the section 5 we focus on the low tem- 
perature regime of the effective 2D model for which a 
RSRG procedure has been developed. 

2. Microscopic model of coupled frustrated spin- 
^ chains 

2.1 Generic Hamiltonian for weakly coupled spin- 
Peierls chains 
We start with the microscopic Hamiltonian 

M L 
a— 1 i—l 

+hi, a Sl a ], (1) 

as a generic model for a quasi one-dimensional frustrated 
spin system 10 with J and a both positive, and S — ^. i 
is a lattice index along the chain of size L and a labels 
the M chains. The effective magnetic field h^ a results 
from a mean-field (MF) treatment 13 of a weak interchain 
coupling, J±Si^ a ■ Si.a+i with Ji <C J, computed in a 
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Fig. 1. Schematic picture of a system of coupled frustrated spin- j 
chains governed by Eq. (1). 



self-consistent way 

h i ,a = J±((Sl a+1 ) + (Sl a _ 1 )). (2) 

The modulation of the nearest neighbor exchange Si^ a is 
also computed self-consistently using 

0~i,a = -^{(Si, a +1 ■ <^i+l,a+l) + (Si. a-1 ' Si+l,a-l))- (3) 
•Jl 

This modulation term might have in fact multiple origins: 
although a four-spin cyclic exchange mechanism provides 
the most straightforward derivation of it, 10 at a qualita- 
tive level, J\ can also mimic higher-order effects in J± 14 
or the coupling to adiabatic phonons. 9 Indeed, in that 
case, due to a magnetoelastic coupling, the modulations 
5i ia result from small displacements of the ions. The elas- 
tic energy is the sum of a local term ^Kn Yl% a ^ia an< ^ 
an interchain contribution K± a 8i >a 8i >a +i of electro- 
static origin. 15 Then, the mean-field equation (3) is re- 
placed by 8 

K\\8%,a + K±(5i,a+1 + £»,a-l) = Jl{Si.a ' Si+l,a), (4) 

giving very similar results 16 so that we can restrict our- 
selves to Eq. (3) without loss of generality. We first con- 
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Fig. 2. SP-AF phase diagram in the (a,Jj_) plane from ED of 
chains of length up to 16 sites. Symbols correspond to different 
values of J4 > as indicated on plot. Typically, finite size effects 
(FSE) arc smaller than the size of the symbols. The computed 
transition lines are extended by tentative transition lines (dashed 
lines) in the region where FSE become large. At J4 = we have 
plotted a few points in the vicinity of the MG point for L = 12 
and L = 16. (Figure reprinted from Ref. 10 ). 



sider the clean undoped system where the M chains 
becomes equivalent (see Fig. 1). By solving the self- 
consistency conditions Eqs. (2-3) using ED on finite 
chains with up to 16 sites 17 (supplemented by a finite 
size scaling analysis), the transition lines Jj_ = J±(a, J4) 
separating the dimerized SP phase (hi^ a = 0) and the 
AF ordered phase (hi, a 7^ 0) have been obtained. 10, 18 
The phase diagram is shown in Fig. 2. If the chains are 
not dimerized, 5Ji^ a is constant and the exchange is just 
renormalised. On the other hand, if the SP phase is sta- 
ble, each chain displays the same dimerized pattern when 
J4 < whereas dimer order is staggered in the transverse 
direction when J4 > 0. Note that, apart from special 
features (see below), physical properties for positive or 
negative J4 are quite similar. Note also that the dimer- 
ized GS would be 2 M -fold degenerate when J4 = (each 
chain is independently two-fold degenerate) while the de- 
generacy is reduced to 2 when J 4 ^ 0. When a = J± = 
but J4 7^ each chain spontaneously dimcrizes and a gap 
opens up. Consequently, a finite value of the AF inter- 
chain coupling is necessary to drive the system into the 
AF ordered phase. 19 The frustration a stabilizes further 
the dimerized phase with respect to the AF one, the crit- 
ical J± (a) increasing with increasing a as seen from the 
phase diagram shown in Fig. 2. 

3. Impurity effects 

Let us now turn to the doped case. A non-magnetic 
dopant is described here as an inert site decoupled from 
its neighbors. Under doping the system becomes non- 
homogeneous so that the MF equations are solved self- 
consistently on finite L x M clusters and lead to a non- 
uniform solution. At each step of the MF iteration pro- 
cedure, we use Lanczos ED techniques to treat exactly 
(although independently) the AI non- equivalent finite 
chains and compute (Sf a ) for the next iteration step un- 
til the convergence is eventually achieved. 

3.1 One impurity: soliton confinement 

We first consider the case of a single dopant. Whereas 
in the case J4 = the soliton remains de-confined, a 
very small J4 7^ is sufficient to produce a confining 
string which binds the soliton to the dopant. 10 In order to 
describe quantitatively the confinement mechanism near 
the dopant, we can define a confinement length in the 
chain direction £j| as 

Eii\(Sf)\ 



(5) 



E,KS?>r 

We have calculated it for a 16 x 8 system with a = 0.5 and 
Ji = 0.1, and we show its variation as a function of J4 in 
Fig. 3. Note that £|| ( J4) 7^ ^ {—Ji) and a power law 5 with 
different exponents rj is expected when J 4 — > 0. A fit gives 
r\ - 0.33 if J 4 < and rj - 0.50 for J 4 > (Fig.3). This 
asymmetry can be understood from opposite renormali- 
sations of J\ for different signs of J4. Indeed, if J 4 < 
then 8Ji : a > and the nearest neighbor MF exchange 
becomes larger than the bare one. Opposite effects are 
induced by J4 > 0. 
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for such a configuration. The behavior of the pairwise 
interaction of two dopants located on different chains 
(Act = 1,2,3) is shown in Fig. 5 for Aa = 1,2,3 for 
J4 > 0. When dopants are on opposite sub-lattices the ef- 
fective interaction is antifcrromagnctic. At small dopant 
separation J eff (Ai) increases with the dopant separation 
as the overlap between the two AF clouds increases until 
Ai ~ 2£||. For larger separation, J eS (Ai) decays rapidly. 
If dopants are on the same sub-lattice, solitons are lo- 
cated on the same side of the dopants and the effective 
exchange J eS (Ai) is ferromagnetic and decays rapidly to 
become negligible when Ai > 2£. The key feature here 
is the fact that the effective pairwise interaction is not 
frustrating (because of its sign alternation with distance) 
although the frustration is present in the microscopic un- 
derlying model. Our next step is to fit the numerical data 



Fig. 3. ED data of the soliton average position vs J4 calculated 
for a = 0.5 and J± = 0.1. Different symbols are used for L X M 
= 12 X 6 and 16 X 8 clusters. The long-dashed line is a power-law 
fit (see text). Inset shows the magnetization profile in the doped 
(a = 1) chain at J4 = 0.08, ie £y ~ 2.5 (Figure reprinted from 
Rcf. 10 ). 



3.2 Two impurities: effective interaction 

We now turn to the investigation of the effective in- 
teraction between dopants. Each impurity releases an ef- 
fective spin i, localized at a distance ~ £y from it due 
to the confining potcnticl set by J4. When two impuri- 
ties are introduced in the system of coupled chains (see 
Fig. 4), we can define an effective pairwise interaction 
J cS as the energy difference of the 5=1 and the 5 = 
ground state (GS). When J cS = E(S = 1) - E(S = 0) is 
positive (negative) the spin interaction is AF (ferromag- 
netic). Let us first consider the case of two dopants in 




Fig. 5. Magnitude of the effective magnetic coupling between two 
impurities located either on the same chain (a) or on different 
ones (b-c-d) vs the dopant separation Ai in a system of size 
L x M = 16 x 8 with a = 0.5, J ± = 0.1, and J 4 = 0.08. Closed 
(resp. open) symbols correspond to AF (F) interactions. Full 
lines are fits (see text). (Figure reprinted from Ref. 18 ). 
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Fig. 4. Schematic picture of the system of coupled frustrated 
spin- i chains doped with two non magnetic impurities separated 
by a distance Ai (resp. Aa) in the chains direction (resp. trans- 
verse direction). 



the same chain, (i) When the two vacancies are on the 
same sub-lattice the moments experience a very small 
ferromagnetic J cff < as seen in Fig. 5 with Aa = 
so that the two effective spins | are almost free, (ii) 
When the two vacancies sit on different sub-lattices, Ai 
is odd and the effective coupling is AF with a magnitude 
close to the singlet-triplet gap. Fig. 5 with Aa = shows 
that the decay of J eS with distance is in fact very slow 



in order to derive an analytic expression for J e g. Using 
only five parameters, two energy scales and three length 
scales, we can fit ED data with very simple mathemati- 
cal expressions. When Aa = (same chain), J cff approx- 
imately fulfills J cff (Ai,0) = J (l - Ai/£j?) for Ai even 
and Ai < (§, and J cff (Ai,0) = other wise. For dopants 
located on different chains and on the same sub-lattice 
(Ai + Aa even) one has, 

J cff (Az,Aa) = - j; c xp(-^)cxp(-^), (6) 

while if the dopants are on opposite sub-lattices, one gets 

J ee (Ai,Aa) = J^exp(-^) (7) 

for Ai < 2£|| and 

J off (Az, Aa) = -f exp(- A '~ 2e|1 ) cxp(-^), (8) 

qi u. 

for Ai > 2£||. The fitting parameters are Jq = .52, 
Jo = 0.3, £{? = 17.33, f || = 2 .5 and ^ = 1 in the case 
considered here : a = 0.5, J j_ = 0.1 and J 4 = 0.08 (see 
Fig.5). 
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4. Impurity induced antiferromagnetic ordering 

In order to study the impurity induced AF ordering 
in the system of weakly coupled chains [Eq. (1)], a long- 
range non-frustrated effective model of diluted effective 
spin- 2 



1 is defined, 



n 



eff 



J2 Cr^^ri-raJS^-S, 



(9) 



with e r = 1 (0) with probability x (1 — x), where x is the 
dopant concentration. Such a model can be studied by 
Quantum Monte Carlo (QMC) simulations on L x x L y 
clusters much larger than those accessible to ED 9 and at 
all temperatures. 
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Fig. 7. Staggered magnetization per site. Disorder average has 
been done over at least 2000 samples, (a) Finite size extrapo- 
lations (see text) at fixed doping x. The different symbols arc 
used for different concentrations x with, from top to bottom, 
x = 4/49, 0.0625, 4/81, 0.04, 4/121, 1/36. (b) Doping depen- 
dance of jtiaf for various numbers of spins and in the thermo- 
dynamic limit (full symbols). (Figure reprinted from Ref. 11 ) 
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Fig. 6. Schematic picture of a doped SP system. Thick bonds 
correspond to dimers, and non-magnetic dopants (released spin- 
i) are represented by open circles (arrows). 



We study the effective Heisenberg model using the 
powerfull Stochastic Series Expansion (SSE) method 20 
to investigate GS as well as finite T properties. In this ap- 
proach, the interactions are sampled stochastically, and 
for a long-ranged interaction the computational effort is 
then reduced from ~ iV 2 to N s In (N s ). 21 In order to ac- 
celerate the convergence of the simulations at the very 
low temperatures needed to study the ground state, we 
use a /3-doubling scheme 22 where the inverse tempera- 
ture is successively increased by a factor 2. Comparing 
results at several f3 = 2™, one can subsequently check 
that the T — > limit has been reached. 

4-1 T = antiferromagnetic ordering 

The AF ordering instability is signalled by the diver- 
gence with system size of the staggered structure factor, 

s(^) = -^—{(j2(-iys!) 2 ). (io) 

L/ X L/y 

Note that within our effective model approach, only 
the sites carrying a "dopant spin" contribute to this 
sum. It is convenient to normalized S with respect to 
the number of sites, i.e. to define a staggered structure 
factor per site; s(ir,ir) = S(n, ir)/L x L y . In an ordered 
AF state, s(ir, 71") should converge, with increasing size, 
to a non-zero value < 1/4. The (finite size) sublatticc 
magnetization ttt-af can then be obtained by averaging 
s(-7r, 7r) over a large number of dopant distributions, i.e. 
(toaf) 2 = 3(s(7r, 7r))di s , where the factor 3 comes from 
the spin-rotational invariance 23 and (. . .)di S stands for 
the disorder average. The staggered magnetization per 
dopant is then simply m sp j n = tuaf/x. Since, strictly 



speaking, in 2D the divergence of S(it, n) occurs only 
at T = (/3 = oo ) it is appropriate to first extrapo- 
late the finite size numerical data to T = 0, and then 
using a polynomial fit in 1/ \fN s an accurate extrapola- 
tion to the thermodynamic limit is performed as shown 
in Fig. 7(a). The doping dependance of the extrapolated 
vtiaf is given in Fig. 7(b). We have tested various fits to 
the data. Assuming a power law oc x^ , the best fit (solid 
line in Fig. 7(b)) gives an exponent \x ~ 1.38 > 1. 

4-2 Neel temperature 

It is also very instructive to calculate the Neel 
temperature, asssuming a small (effective) 3D magnetic 
coupling A30 between the 2D planes. Using an RPA cri- 
terion, the critical temperature Tjv is simply given by 
Xstag(Tjv) = 1/1 where the staggered spin suscepti- 
bility (normalized per site) is defined as usual by, 

X, t jT) = l^J2(- 1 ) n+rj f dT{S*(0)S?(r)) , (11) 

and averaged over several disorder configurations (typi- 
cally 2000). Since Xsta g (T/v) is expected to reach its ther- 
modynamic limit for a finite linear size L as long as T/v 
remains finite, accurate values of Tjv can be obtained us- 
ing a finite size computation of Xstag(T") for not too small 
inter-chain couplings. Fig. 8(a) shows that Xstag^) di- 
verges when T — > 0. X\r is determined by the intersec- 
tion of the curve Xstag (T) with an horizontal line at co- 
ordinate I/A3D. Note that finite size corrections remain 
small, even in the worst case corresponding to very small 
Xzd values and large dopant concentrations. The doping 
dependance of T/v is plotted in Fig. 8(b) for a particularly 
small value A3D = 0.01 (in order to show the small size 
dependance observable in that case). It clearly reveals a 
rapid decrease of T/v when x — > 0, but, in agreement with 
experiments, does not suggest a non-zero critical concen- 
tration. In Fig. 8(b), we show the behavior of T/v(x) down 
to x ~ 0.007. Note that from numerical fits of the data, 
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Fig. 8. (a) Staggered susceptibility of a 2D layer vs temperature 
(using log-log scales) for N s = 256 (full symbols) and N s = 576 
(open symbols) spins. Concentrations x are shown on the plot. 

= 100 is shown by the dashed line, (b) Neel temperature 
vs dopant concentration x for a 3D RPA inter-plane coupling 
\ ZD = 0.01 and for N s = 256 and N s = 576 spins. (Figure 
reprinted from Rcf. 11 ) 



we can not clearly distinguish between a power-law be- 
havior (with an exponent ~ 2.5) and an exponential law 
like Aexp(— B/x). as suggested by fits of experimental 
data for Cui_ x Zn K Gc0 3 . 24 

5. Low temperature critical properties 

5. 1 Saturation of the Curie constant 



As a first investigation, we have computed the uni- 
form susceptibility x(T) for a wide range of temperatures 
using SSE simulations. Results for the Curie constant 
C(T) = Tx(T) are shown in Fig. 9. At the highest tem- 
peratures the effective dopant spins behave as free spins 
while at low temperature we observe a new Curie- like be- 
havior with a reduced Curie constant ~ 1/12. Although 
here the 2D system orders at T = (as proven above) 
this behavior agrees with a qualitative argument based 
on the formation of large spin clusters due to the presence 
of AF as well as ferromagnetic (F) couplings, applied by 
Sigrist and Furusaki 25 in the ID case. 

5.2 Real Space Renormalization Group 

In order to explore this large spin formation, we fol- 
low the pioneering work of Ma, Dasupta and Hu 26 and 
also Bhatt and Lee, 27 and extend the Real Space Renor- 
malization Group (RSRG) scheme to hamiltonian (9) 
with F and AF long-distance couplings. It turns out that 
such a RSRG procedure has allowed the identification 
of the so-called Infinite Randomness Fixed point (IRFP) 
by D. Fisher for the disordered Heisenberg spin-i chain 28 
as well as for the Ising chain in transverse field at crit- 
icality. 29 Nevertheless, such an IRFP (as well as its as- 
sociated random Singlet Phase) is unstable against the 
introduction of randomness in the couplings' signs, as 
shown first by Westerberg et al. 30 for the Heisenberg 
spin-i chain with random couplings in magnitude and 
sign. The new fixed point identified in the ID case, using 
an extended RSRG scheme 31 and checked numerically 
using QMC simulations, 32 is especially characterized by 
a large spin formation. 33 
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Fig. 9. Curie constant per spin (C) plotted vs the energy scale. 
Quantum Monte Carlo SSE results shown vs T for N s = 256 
spins and 9 different concentrations x indicated on the plot. Dis- 
order is performed over 10 3 to 10 4 samples. (Figure reprinted 
from Rcf. 12 ). 



Fig. 10. Schematic picture of one RG step performed in the 
doped spin-Peierls system. Thick bonds stand for dimers and 
non-magnetic impurities (released spins i) are represented by 
open circles (black arrows). The initial RSRG step is illustrated 
starting from a typical configuration with 4 impurities: (a) The 
strongest coupled pair is identified (red line), (b) This pair, e.g. 
ferromagnetic here, is replaced by a spin 5 = 1 (red arrow), (c) 
The couplings with all other spins (dashed lines) are renormal- 
ized. 



In the 2D problem adressed here, long range inter- 
actions with random sign and magnitude can also be 
studied by an extended version of the RSRG technique. 
Let us first define the effective interaction as J;.j where 
i and j label the randomly distributed spins and run 
from 1 to N s . One single RG step is described as fol- 
lows (see Fig. 10): 1) Identify the most strongly coupled 
pair of spins (SijSa) i-c. with the largest energy gap 
Ai, 2 , Ai, a = Ji, 2 (l + \Si -S 2 |) if Ji,2 > (AF) and 
Ai, 2 = -Ji, 2 (5i + S 2 ) if Ji,2 < (F). Note that A li2 
defines the energy scale of the transformation which will, 
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in fact, play the role of the temperature. 2) Replace it by 
an effective spin S = |Si — S 2 | if the coupling is AF or 
S = Si + S 2 in the F case. 3) Renormalize all the mag- 
netic couplings with the following rules : (i) If S 7^ 0, as 
given by a first order perturbation theory, the new cou- 
plings between S and all the other spins (S3, S4, Sn s ) 
are set to 



with 



(S',S t ) 



Ji,i c(Si, S 2 ,S ) + J 2 .i c(S 2 ,S 1 ,S ), (12) 



S'{S' + 1) + SijSi + 1) - S 2 (S 2 + 1) 
25" (£" + 1) 



c(Si, S 2 , 5 ) 

(13) 

(ii) If S =0, the pair (Si,S 2 ) is frozen. Using a clus- 
ter approximation 27 that involves only the extra pair 
(53,54) the most strongly coupled to Si and S 2 and a 
second order perturbation, the coupling ^4 is renormal- 
ized as 



3,4 



J3A + 



2Si(Si + 1) 
3Ji.2 



(Jl,3- J2, 3 )(J 2 ,4- Jl.A). (14) 



The same procedure is then reiterated. We also check 
that the RSRG preserves the non-frustrated character of 
the problem. 

5. 3 Large spin formation 

Due to the presence of both F and AF couplings, clus- 
ters with large effective spins are created during the pro- 
cedure similarly to what occurs in the ID random F- 
AF spin-i chain. 30 At each RG step, the energy scale 
Ao decreases and both the number of inactive spins 
frozen into singlets and the number of clusters build 
from a large number n correlated spins- \, increase. The 
aforementioned random walk argument predicts that, 
the average number (n) of spins- h inside clusters and 
their average spin magnitude (S ) should be related by 
(5 eff ) ~ (n) 1 / 2 at low enough temperatures. Therefore 
we expect the effective spin of these clusters to grow 
monotonously as the energy scales down. We have an- 
alyzed this process using the RSRG scheme to compute 
both (S eS ) and (n) as a function of (Ao). This is shown 
in Fig. 11 which clearly demonstrates the formation of 
large moments. Moreover, power-law divergences like 

(5 cff ) ~ (A )- a ^ and (n) ~ (A y K(x) (15) 

are observed with k ~ 2a. Interestingly enough, we find 
that a depends on x, like ~ y/x, in contrast to the ran- 
dom F-AF spin chain for which a = 0.22 ± 0.01. 7 ' 31 ' 32 

The Curie constant can also be evaluated using the 
decimation procedure. The RSRG computation of C is 
performed, at each RG step, using the formula C = 
zjj-^Za N a o-(<7 + 1) where N a is the number of active 
effective spins of size a, the data being then averaged 
over disorder. We have also checked that finite size ef- 
fects are negligible when N s > 256 and we have chosen 
N s = 1024 in most computations. Fig. 12 shows results 
for the Curie constant per spin obtained with the deci- 
mation procedure. We can confront the RSRG results to 
QMC data obtained on the same model and for the same 
parameters. We observe a qualitative agreement between 




Fig. 11. Average effective spin (S eff ) of the clusters of active spins 
vs the energy scale (Ao) for six different concentrations x indi- 
cated on the plot. Numerical RSRG data obtained for N s = 1024 
spins over more than 10 4 samples. Inset: for the same samples, 
average number (n) of initial spins- i per cluster vs (Ao). Dashed 
lines are power-law fits (see text). (Figure reprinted from Ref. 12 ). 
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Fig. 12. Curie constant per spin (C) plotted vs the energy scale. 
Numerical RSRG results shown for N s = 1024 spins and 6 differ- 
ent concentrations x as indicated on the plot, vs the RG energy 
scale (Ao). Error bars are typically smaller than symbol sizes, 
the number of samples always exceeding 10 4 . The full line cor- 
respond to the saturation value of 1/12. (Figure reprinted from 
Ref. 12 ). 



Fig. 9 and Fig. 12 where, at high temperatures, the spins 
behave as paramagnetic free magnetic moments (giving 
a Curie constant of j per spin) and where saturation to 
J2 is observed at low T, the spins being correlated inside 
large clusters. At small concentration, the agreement be- 
comes even quantitative, as seen in the Fig. 13 where a 
direct comparison between RSRG and SSE is shown for 
the six lowests values of x. 
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Fig. 13. Direct comparisons between RSRG (full lines) and SSE 
simulations (symbols) of (C) — 1/12 are shown for the 6 different 
concentrations indicated by (a), (b),..., (f) vsthe RG energy scale 
(Aq) or the SSE temperature T. (Figure reprinted from Ref. 12 ). 



5.4 Scaling regime 

We now turn to the analysis of the scaling regime. At 
very low temperatures, from the analogy with ID sys- 
tems we expect the following quantum corrections 

(C(T)) — 1/12 ~ T 7 , (16) 

with 7 — a in ID. 32 Similarly, one expects for the stag- 
gered susceptibility (per spin), 



(17) 



In ID, 7' is expected to be equal to a, 35 but in the present 
case, by direct fits of the low T (see Fig. 14), we found 

7' ~ 7 ~ 2a. 

Thoses exponents are plotted in Fig. 15. An overall very 
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Fig. 14. T(xstag(T)) plotted vs T for five different concentra- 
tions. Full lines are fits corresponding to power-law behaviors 
~ T~ 2 ~t' . All data are computed by QMC using N„ = 256 ran- 
dom spins and averaged over disorder. (Figure reprinted from 
Rcf. 12 ). 



good agreement is seen between the different methods, 
in particular between the estimates of 7 obtained from 



the analysis of the Curie constant computed by QMC 
and RSRG. We stress again that the exponent a de- 
duced from the analysis of the change of cluster sizes 
and spins is roughly a factor of 2 smaller than 7. Inter- 
estingly enough, 7 ~ 2a cx \fx in all cases. 
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Fig. 15. Exponents a(x), "y(x) and 7' (2;) extracted from various 
SSE and RSRG data as indicated on the plot (see text for de- 
tails). Straight lines arc ~ ^Jx. (Figure reprinted from Ref. 12 ). 



It is remarkable that such properties are observed 
above the T = ordered magnetic groundstate. More- 
over, it is the first example of a two- dimensionnal ran- 
dom magnet exhibiting a large spin phase with a dis- 
order (the concentration x) dependence of the critical 
exponents. Whereas doped CuGeC>3 might be a good 
candidate for the experimental observation of a critical 
regime, the so large three-dimensionnal ordering temper- 
ature prevents such an observation. More strongly di- 
luted samples are necessary to reach the scaling regime. 

6. Conclusion 

Starting from a microscopic model for impurity-doped 
dimerized and frustrated coupled spin-chains (similar to 
let's say, Zn-dopcd copper gcrmanate), we have derived 
a simple and tractable effective model describing S=l/2 
effective spins interacting through long range (spatially 
anisotropic) interactions. An important property of this 
long-range exchange is its oscillating nature (with both 
AF and F couplings) in such a way it is essentially non- 
frustrating, hence leading to an AF instability at low 
temperature (in the presence of a small 3D coupling). 
This is a remarkable fact since magnetic ordering in the 
pure system is prevented by frutration. Using state-of- 
the-art SSE QMC computations thermodynamic prop- 
erties of the effective model are studied in details down 
to very low temperatures. Critical behaviors with impu- 
rity concentration-dependent exponents are found and 
interpreted using a new Real Space RG procedure clearly 
establishing the important role of large clusters of corre- 
lated spins carrying a finite magnetization. 
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